LAPTH-044/11 



Future constraints on neutrino isocurvature perturbations in the curvaton scenario 

Eleonora Di Valentino/ Massimiliano Lattanzi,-^ Gianpiero Mangano,'^ Alessandro Melchiorri,^ and Pasquale Serpico'* 

^Physics Department and INFN, Universitd di Roma "La Sapienza" , Pie. Aldo Mora 2, 00185, Rome, Italy 
^ Dipartimento di Fisica G. Occhialini, Universitd, Milano-Bicocca and INFN, 
sezione di Milano-Bicocca, Piazza della Scienza 3, 1-20126 Milano, Italy 
^ INFN, Sezione di Napoli, Complesso Univ. Monte S. Angelo, Via Cintia, 1-80126 Napoli, Italy 
'^LAPTh, Univ. de Savoie, CNRS, B.P.llO, Annecy-le-Vieux F-74941, France 

(Dated: March 13, 2012) 

In the curvaton scenario, residual isocurvature perturbations can be imprinted in the cosmic 
neutrino component after the decay of the curvaton field, implying in turn a non-zero chemical 
potential in the neutrino distribution. We study the constraints that future experiments like Planck, 
SPIDER or CMBPol will be able to put on the amplitude of isocurvature perturbations in the 
neutrino component. We express our results in terms of the square root 7 of the non-adiabaticity 
parameter a and of the extra relativistic degrees of freedom AA'cff • Assuming a fiducial model with 
purely adiabatic fluctuations, we find that Planck (SPIDER) will be able to put the following upper 
hmits at the la level: 7 < 5.3 x 10"^ (1.2 x 10"^) and AA^eff < 0.16(0.40). CMBPol will further 
improve these constraints to 7 < 1.5 x 10~^ and AAeff < 0.043. Finally, we recast these bounds 
in terms of the background neutrino degeneracy parameter f and the corresponding perturbation 
amplitude a^, and compare with the bounds on ^ that can be derived from Big Bang Nucleosynthesis. 



PACS numbers: 98.80.Cq, 98.70.Vc, 98.80.Es 

I. INTRODUCTION 

In single-field inflationary models, the same field is re- 
sponsible for driving an accelerated expansion stage and 
for the generation of a nearly scale invariant primordial 
perturbation spectrum. As there is only one degree of 
freedom, this class of models predicts that perturbations 
are necessarily adiabatic, i.e., the ratio between the num- 
ber densities of the different particle species is spatially 
homogeneous. Significant non-Gaussianities in the fluc- 
tuations are also excluded. To date, both these features, 
adiabaticity and gaussianity, are consistent with data. 

However, the presence of a significant, albeit sub- 
dominant, non-adiabatic (otherwise called isocurvature) 
perturbation component cannot be excluded, see e.g. [1]- 
[Bj. This component must be necessarily related to some 
extra field other than the infiaton, as in multifield in- 
flationary models, where non trivial trajectories in field 
space are possible. Since in these cases the adiabatic 
and isocurvature fiuctuations would be related to differ- 
ent fields, generating a sizeable isocurvature fluctuation 
requires in general a certain amount of flne-tuning. 

A different mechanism for isocurvature modes produc- 
tion has been proposed in [71 [5] . While the inflaton is 
only responsible for driving the exponential expansion, 
primordial fluctuations are generated by a "curvaton" 
fleld. The initial isocurvature perturbation in the cur- 
vaton is then converted into an adiabatic component af- 
ter inflaton decay. This model allows for some residual 
isocurvature components imprinted in the other compo- 
nents of the cosmological fluid, cold dark matter, baryons 
and neutrinos, after curvaton decay. In particular neu- 
trino isocurvature perturbations requires a non vanish- 
ing chemical potential for their background distribution 
in phase space. Probing their non adiabatic perturba- 



tions is thus, a way to constrain the lepton number in 
neutrino sector. An analysis of the bounds on neutrino 
isocurvature perturbations using recent data is presented 
in Ref. while limits on isocurvature perturbations in 
an extra radiation component (not necessarily related to 
neutrinos) have been derived in Ref. [10] . Quite inter- 
estingly, a future detection of isocurvature perturbations 
will allow for a reconstruction of the inflationary poten- 
tial 

The aim of the present paper is to assess the capability 
of future cosmic microwave background (CMB) experi- 
ments like Planck [HI [13], SPIDER [14^ and CMBPol 
[15) to constrain simultaneously the amplitude of isocur- 
vature perturbations in the neutrino component and the 
extra energy density associated to the neutrino chem- 
ical potential. The bounds can then be translated into 
constraints on the neutrino chemical potential to temper- 
ature ratio {i = e, fj,, r) and the corresponding pertur- 
bation amplitudes. These are complementary to bounds 
on the ^i's which can be derived using Big Bang Nucle- 
osynthesis (BBN). Light nuclei yields in fact, are quite 
strongly influenced by neutrino asymmetries, mainly in 
the I'e sector, see e.g. [TB|. Since flavour oscillations 
are efficient in mixing different fiavour distributions, the 
three parameters are driven to almost the same value at 
the onset of BBN [T7] , with possible differences in the e 
and ^, T sectors which cannot be larger than few percents 
depending on the value of the 6113 mixing angle [HI \W\ . 

The paper is organized as follows. In Section [llj we re- 
view the neutrino isocurvature perturbations which are 
generated in the curvaton scenario. Section |III| contains 
a forecast analysis of bounds on these perturbations from 
future experiments, while a comparison with the corre- 
sponding BBN constraints is described in Section |TV| Our 
conclusions are reported in Sec. [Vj 
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II. NEUTRINO ISOCURVATURE 
PERTURBATIONS 

Density perturbations are conveniently described in 
terms of the gauge-invariant quantity ^ 



p 



(1) 



where ij) is the (gauge-dependent) curvature perturba- 
tion, H the Hubble parameter, p the total energy den- 
sity, and the dot denotes derivatives with respect to the 
cosmological time t. 

The quantity C describes the curvature perturbation on 
slices of uniform total density. In the case of multicompo- 
nent fluids, it is useful to define quantities Ci describing 
the curvature perturbation on slices of uniform density 
of the i-th component 
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When dealing with cosmological neutrinos, it is cus- 
tomary to define the effective number of neutrino families 
iVoff as the ratio between the total neutrino density and 
the density of a single non-degenerate = 0) neutrino 
species in thermal equilibrium at Tiy — (4/ll)^/^T^. In 
the standard cosmological scenario iVoff — 3.046, see [24], 
and any deviation ANce from this value indicates the 
presence of an extra energy density of rclativistic parti- 
cles in the early Universe. It is clear, from our definition, 
that N^ff = We can thus relate the isocurvature 

perturbation in the total neutrino density to the fluctu- 



ations 6N^ 



W. 



eff ■ 



Si, — 3(Ciy — C7) 
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An adiabatic fluctuation is defined as one for which 
the ratios Spi/pi are all the same, so that Q = C for all 
components. Correspondingly, a nonadiabatic (or isocur- 
vature) fluctuation Si in the i-th fluid component is de- 
fined as the relative entropy fluctuation with respect to 
photons: 

S, EE 3(C. - C7) ■ (3) 

In the following, we shall consider neutrinos with an 
equilibrium distribution function 

f,iE)^[eME/nT^^)^^ , (4) 

where T^, is their temperature, and — pi/T^, pi being 
the chemical potential. The index i runs over the three 
standard model neutrino families, i — e, p, r, and the 
minus (plus) sign is for neutrinos (antineutrinos). No- 
tice that the existence of neutrino isocurvature pertur- 
bations necessarily implies a non zero lepton asymmetry 
in the neutrino sector, = — np, unless the asymme- 
tries in the three flavours exactly cancel. At this stage, 
we have allowed for the possibility of the three neutrino 
families having different chemical potentials. The neu- 
trino temperature is T^, — until the time of electron- 
positron annihilation, occurring at ~ 1 MeV (shortly 
after neutrino decoupling), while at later times it is given 
by Ti, = (4/ll)^/'^T-y, up to tiny corrections due to neu- 
trino reheating at the annihilation stage [24,. 

Given the distribution function Eq. Q, the energy 
density pi = p^. + pp. in the high-tcmperature limit 3> 
iTii, writes [25] : 



III. CMB CONSTRAINTS AND FORECAST 

In the following, lacking a better theoretical motiva- 
tion, for simplicity we shall assume that both the average 
values and the fluctuations in the chemical potentials are 
flavor blind, i.e. £,& = = £,t = and similarly for the 
(5^'s. Also, we assume that fluctuations in the neutrino 
degeneracy parameter are gaussian distributed with vari- 
ance cr| around the mean ^. In general, both quantities 
can have a scale and epoch dependence. 

Conventionally, rather than in terms of Si, of Eq. ([7| , 
in CMB studies the "non-adiabaticity" of perturbations 
is expressed in terms of the ratio of the power spec- 
trum Ps{k) of isocurvature perturbations to the curva- 
ture perturbation spectrum P^(fc), evaluated at a fixed 
pivot wave number fco — 0.002 Mpc^^. In particular, one 
introduces the quantity a defined by [U [5] 

«(fco) ^ Psikg) , . 

1 - a{ko) ~ Pciko) ' ^' 

Another necessary ingredient to be taken into account 
is the correlation between the adiabatic and isocurva- 
ture modes p5H28j . Given the cross-correlation power 
spectrum P^5(A:), this is parameterized in terms of the 
cross-correlation coefficient /3, defined as 

B = (9) 

VPsiko)Pdko) 

We remark that we choose the sign convention for the 
curvature perturbation such that the temperature fluctu- 
ation at large scales is given by AT/T = (/5 — 25/5. In 
terms of the variables used in the WMAP analysis [61129]. 
(■ = ??. = —7?,, and our definition of /3 coincides with the 
one used there. In this case, the physically observable 
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effect is that correlated perturbations {f3 > 0) reduce the 
temperature power spectrum at low multipoles. 

Given the above convention, the adiabatic and isocur- 
vature fluctuations in the curvaton scenario are totally 
anticorrelated so that in the following we will 

always take /3 = — 1. We also take the two power spectra 
to have the the same spectral tilt rig: 



27r2 



(10) 



Note that the CMB is sensitive to the parameters of 
the scenario not only via a, but also via the total AiVofj 
induced by the average value ^ and, if sufficiently large, 
in principle also by the variance (t|. For analyses or 
forecasts, one should thus consider the constraints in the 
A-^eff — 01 plane. 

Current WMAP7 bounds on totally anticorrelated 
isocurvature perturbations are at the level a(A;o) < 1.1 x 
10~^ (at the 95% confidence level); inclusion of addi- 
tional datasets can improve this bound by a factor 2 
or 3 [21 [33]. For what concerns the effective number 
of relativistic species, WMAP7 observations only pro- 
vide a lower limit iVoff > 2.7; interestingly, when other 
cosmological measurements are considered, the result is 
A'cff = 4.34 ± 0.9 [6. , indicating that the data seem to 
prefer ANes > 0. 

The total CMB power spectrum can be parameterized 
in terms of the adiabatic, neutrino isocurvature density 
and totally anticorrelated spectra as follows 



with a defined in Eq. 



(11) 

A shortcoming of this pa- 
rameterization is that the partial derivative needed 
for the Fisher matrix computation (see below), diverges 
for a = 0. This prevents the use of the Fisher matrix 
formalism for the fiducial value a = 0. For this reason, 
we find convenient to introduce the auxiliary parameter 
7 — \foL and write 



-27v/(i - -r^)cr" ■ 



(12) 



It is straightforward to check that the partial derivative 
^ is finite for 7 = 0. 

In the following we will derive forecasts for the Planck 
[12, SPIDER [H] and the CMBPol [H] experiments. 
The Planck satellite [H] US], launched in May 2009, is 
currently measuring the CMB temperature and polariza- 
tion fluctuations with unprecedented precision {AT/T ~ 
2 X 10~^) over the whole sky and down to very small 
angular scales (~5'). Planck measurements, planned to 
be publicly released to the scientiflc community in Jan- 
uary 2013, will signiflcantly improve the determination 
of cosmological parameters and will allow to test further 
the ACDM paradigm. SPIDER [M], scheduled to flight 
in 2012, is a ballon-borne polarimeter design to accu- 
rately measure the i?-mode of CMB polarization down 



to £ ~ 100. Finally, CMBPol [TS] is a next-generation 
satellite currently in the concept study phase. 

In order to derive forecasts for these experiments, we 
use a Fisher matrix formalism, for three frequency chan- 
nels for each experiment (the experimental specifications 
are listed in Table |l]). We consider a detector noise of 



Experiment Channel [GHi 


^] FWHM 






Planck 


217 


5.0' 


13.1 


18.5 


fsky = 0.65 


143 


7.0' 


5.99 


8.48 




100 


9.5' 


6.75 


9.55 


SPIDER 


280 


17' 


0.20 


0.29 


fsky =0.1 


150 


30' 


0.08 


0.11 




90 


49' 


0.08 


0.11 


CMBPol 


220 


3.8' 


0.66 


0.93 


fsky = 0.65 


150 


5.6' 


0.25 


0.35 




100 


8.4' 


0.22 


0.31 



TABLE I. Experimental specifications for Planck [T5], SPI- 
DER [13 and CMBPol [15]. For each experiment, we list the 
observed fraction fsky of the sky, the channel frequency in 
GHz, the FWHM in arcminutes, the sensitivity per pixel for 
the Stokes / (or), Q and U (crp) parameters in ^K. 



(^cr)^ for each frequency channel where d is the FWHM 
of the beam assuming a Gaussian profile and a is the sen- 
sitivity. We add to each fiducial spectrum Ce, calculated 
with CAMB [51], a noise spectrum given by 



(13) 



where lb = V& In 2/6*. In the analysis, we assume that 
beam and foreground uncertainties are smaller than the 
statistical errors. 

The Fisher matrix is defined as 



(14) 



where £(data|p) is the likelihood function of a set of pa- 
rameters p given some data; the partial derivatives and 
the averaging are evaluated using the fiducial values po 
of the parameters. The Cramer-Rao inequality implies 
that {F^^)ii is the smallest variance in the parameter 
Pi, so we can generally think of as the best possible 
covariance matrix for estimates of the vector p. The 1-a 
error for each parameter is then 



<^P. = ViF-^h ■ (15) 
The Fisher matrix for a CMB experiment is (see [35]) 



CMB 



2J 



1=2 



dpi ' 



(16) 



where Cov; is the spectra covariance matrix. We use 
information in the power spectra up to Zmax — 2500. The 
partial derivative is analytical in 7 = 0: 



d-f 



(-27)Q° 



ad 



2(1 - 27^) 

y(w) 



(17) 
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As anticipated above, the parameterization in terms of 7, IV. COMPARISON WITH BBN CONSTRAINTS 

instead of a = 7^ as often seen in the hterature, cancels 

Big Bang nucleosynthesis, and in particular the pri- 
mordial helium abundance Yp, is recognized to be the 
most sensitive cosmic "leptometer" presently available, 
see for example [3S1 [37] or the review [TB] . So, it would 
be interesting to compare BBN constraints to the ones 
derived above. This task is made non-trivial by the fact 
that BBN is sensitive to different parameters than the 
CMB in particular, to a combination of the role of N^s, 
entering the expansion rate of the universe, and in prin- 
ciple to all the parameters describing the distribution of 
the I'e-flavour neutrinos. In the case of interest, which as- 
sumes flavour-independent parameters and gaussian dis- 
tributions, the only two independent parameters turn to 
be ^ and cr^ , with A A^cff fully specified in terms of them, 
but subleading and essentially negligible for the values 
^ <C 1 of interest here. Even assuming that the av- 
erage value f is scale-independent, a slight dependence 
on the scale is expected for the width of the distribu- 
tion of fluctuations. Let us fix (arbitrarily) at a scale 
Abbn, roughly corresponding to the horizon size at the 
time of BBN, of the order of ^ 0(100) comoving par- 
sees. Namely, we fix (t| ^ A|(A:bbn) where /cbbn = 
Stt/Abbn = 6 X lO^'Mpc^^ The CMB constraints can 
be translated into cr| by just evaluating A^(/cbbn) (given 
that A^ has the same scale-dependence). Using WMAP7 
best fit values A2(fc = 0.002 Mpc~^) = 2.42 x 10"^ and 

= 0.966 gives A^(A:bbn) = 1-35 x 10"^ A first im- 
portant consequence of this estimate is that the order 
of magnitude of the present constraints from CMB on a 
also holds for BBN-relevant fluctuations. In turn, it can 
be seen that this implies that cr^ is very small. This is 
an important information, since it allows us to use the 
predictions of homogeneous, degenerate BBN to infer the 
results of an otherwise inhomogeneous degenerate BBN 
scenario (see [38] for an early study of this subject). In 
fact, for a gaussian probability distribution for ^, 



the divergence of the partial derivative in a = 0. 
Thus this parameterization allows us to use the Fisher 
matrix formalism for the fiducial value 7 = 0. 

In the present analysis, we take as a fiducial model 
a flat ACDM model with parameter values given by 
the WMAP7 measurementf|5 i-e. ^bh"^ = 0.02258 
and Vldmh^ — 0.1109, the optical depth to reioniza- 
tion T = 0.088, iJo — 71km/s/Mpc, the spectral index 
Us = 0.963, and the amplitude of the curvature pertur- 
bation A^(A:o) = 2.43 x 10"^. We consider three families 
of massless neutrinos, but we checked that taking mas- 
sive neutrinos with total mass M^, = 0.6 eV neutrinos did 
not affect the results. Finally, we take the fiducial values 
AiVefl = 0, 7 = 0. 

The results of our analysis are shown in Figure [T] 
where we draw the 2-dimensional likelihood in the AA'off- 
7 plane for Planck, SPIDER and CMBPol. The corre- 
sponding 1-cr constraints for 7 and AA^eg are reported in 
Tab. nil 




0.50 

AN„„ 



0.75 



1 .00 



FIG. 1. 68% and 95% c.l. likelihood contours for Planck (solid 
line), SPIDER (dot-dashed line) CMBPol (dashed line). 



P(e) = (27ra|)-V2exp[-(e-OV(2^|)] , (18) 
one can estimate, for a generic nuclide abundance A, 



fiducial value cr(Planck) cr(SPIDER) a(CMBPol) 



7 

AiVeff 



0.0 
0.0 



5.3-10"^ 1.2-10"^ 
0.16 0.40 



1.5- 10" 
0.043 



TABLE II. l-(7 constraints for 7 and AAofi, for the Planck, 
SPIDER and CMBPol experiments. 



^ http : //lambda. gsfc .nasa.gov/product/map/current/params/ 
r lcdm_sz_lens_wmap7 . cf m 



{X) 



p(e)[A(c) + x'm - + o((e - ifm 



0{al). 



(19) 



The vanishing of the integrand linear in ^ depends on the 
fact that P(0 is an even function of ^ — ^. Additionally, if 
the curvature of the function A(^) is relatively small (as 
it happens to be, see Fig. 13 in [T^) the approximation 
(A) A(f) works even better (see also figs, in [38]). We 
estimated that even for a value as large as ~ 0.1 the 
error of the approximation with respect to a proper av- 
eraging is of ^ 0.6% for deuterium (hence well below the 
observational error) or of the order of 0.3% for helium- 
4, comparable with the theoretical error and well below 
the error on the observations. For smaller cr^, it scales 
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as cr| and becomes soon negligible. As a consequence, 
the bounds computed in homogeneous, degenerate BBN 
can be used, to an excellent approximation, also for the 
case at hand. Needless to say, this also implies that BBN 
may give excellent constraints on i^, but it is insensitive 
to physically relevant values of the fluctuation a^. 

By using the same conservative input as in 39 (fourth 
line in their Table I), we obtain the bounds 



-4.0 



^™" = -0.055, ^'"^'^ = 0.12, 



(20) 



corresponding to the value below/above which only 5% 
of the area of the marginalized distribution of probability 
of ^ lies, respectively. The BBN computation is based on 
the PArthENoPE code gO]. 

In order to illustrate the synergy between BBN and 
CMB, it turns useful to translate the CMB forecasts in 
the ^ ~ (T(^ plane. For simplicity, let us write down the 
relation between variables in the (plausible) assumptions 
of ^ 1 and 6^ <^ Then, the relation between the 
power spectrum of the isocurvature perturbation S,^ to 
that of ^ writes 



TT X 

15 



2?! 



(21) 



The relation above implies A| = F(^)^A|, so that 

AK^bbn) - Fi^fal ^ a ~ 7.4 X 10«F(OV| , (22) 

where we have used the fact that the data constrain a to 
be 0(0.01) or less. Recalling that 7 = ^/a, we finally get 
the relation that we were looking for 



7 ~ 2.7 X 10'*F(^)cr^. 



(23) 



On the other hand, recalling that iVcff — J^^ij with Ai 
given by Eq. ^, we can translate the bounds from the 
(7, AA'eff) plane to the (^, cr^) plane. For this analysis, 
we only consider Planck and CMBPol since they give the 
better constraints on the parameters and can possibly 
become competitive with BBN in this respect. The two- 
dimensional 68% and 95% confidence regions for Planck 
and CMBPol are shown in Fig. [2j alon g w ith the present 
BBN constraints on ^ reported in Eq. ( 20 ) . 



V. CONCLUSIONS 

We have considered models in which a residual com- 
ponent of isocurvature fluctuations, and consequently a 
non-zero chemical potential, are generated in the neu- 
trino component after curvaton decay. Using Fisher ma- 
trix techniques, we have assessed the constraints that 
Planck, SPIDER and CMBPol wih be able to put on 
the amplitude of the isocurvature component and on the 
extra energy density associated to the non- vanishing neu- 
trino chemical potential. These bounds have been ex- 
pressed in terms of the average value of the degeneracy 
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FIG. 2. 68% and 95% c.l. likelihood contours in the 
(log j^j , log (Tj) plane for Planck (solid line) and CMBPol 
(dashed line). The BBN allowed region are also shown (left 
of v erti cal lines), corresponding to the case considered in 
Eq. pOl. 



parameter f and its spatial variance cr|, and then com- 
pared with the constraints resulting from the effect of 
neutrino degeneracy on BBN. While the latter is only 
sensitive to the mean value ^, for small ct^, CMB data 
provide a negative correlation between ^ and a^: for 
(5^ ^ ^ (on which our analysis is based), large values 
of the fluctuations are allowed for sufficiently small 

In particular, assuming a flducial model with purely 
adiabatic primordial fluctuations, we find that the future 
experiments will sensibly improve the constraints on the 
non-adiabaticity parameter 7 and on the effective num- 
ber of neutrino families AA'cff. The current 95% C.L. 
WMAP bound corresponds to 7 < 0.1; we find that, 
at the same level, SPIDER will be able to constrain 7 
below 2.4 • 10~^, representing an improvement of a fac- 
tor 4. Planck and CMBPol wiU be able to put 95% 
C.L. upper limits 7 < 1.1 • 10"^ and 7 < 3.0 • 10"^, 
i.e. to improve by a factor 10 and 30 over current data, 
respectively. For what concerns AiVeff, we already no- 
ticed how WMAP only provides a lower bound for this 
quantity. In general, current data allows to constrain 
this parameter with precision cr(AA^off) — 1 or smaller 
(see e.g. Ref. [H] for a detailed analysis). In terms 
of the average value of the degeneracy parameter, this 
read^ a{^) ~ 0.9, which is quite large with respect 
to the BBN bound. For comparison, we find that if 
ATVcff = 0, Planck, Spider and CMBPol wfll be able 
to bound AiVoff < 0.3, 0.8,_0.08 at the 95% C.L., re- 
spectively, corresponding to ^ < 0.5, 0.8, 0.24. Although 



^ We note that for finite neutrino masses the effects of f on the 
CMB observables are not completely encoded by ATVofj [251 142] . 
Joint constraints on ^ and AN^ff have been derived for example 
in Refs. ,42..43| . 
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these values represent considerable improvements over 
present CMB constraints, they show that future CMB 
experiments, with the partial exception of CMBPol, will 
still be unable to compete with BEN in this respect. 
Moreover BBN is the only one sensitive to the sign of ^. 
Both effects are actually due to the dominant weak in- 
teraction probe provided by BBN, as opposed to purely 



gravitational effect to which CMB is sensitive. 
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